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We investigate the parity-broken phase (Aoki phase) for staggered- Wilson fermions by using the 
Gross-Neveu model and the strong-coupling lattice QCD. In the both cases the gap equations 
indicate the parity-broken phase exists and the pion becomes massless on the phase boundaries. 
We also show we can take the chiral and continuum limit in the Gross-Neveu model by tuning 
mass and gauge-coupling parameters. This supports the idea that the staggered- Wilson fermions 
can be applied to the lattice QCD simulation by taking a chiral limit, as with Wilson fermions. 
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1. Introduction 

Recently staggered-based Wilson fermions were proposed by introducing the taste-splitting 
mass or the flavored-mass terms into staggered fermions [III 0, Hj. They can be applied to lattice 
QCD not only as Wilson fermions but also as an overlap kernel. One possible advantage of these 
novel fermions called staggered-Wilson and staggered-overlap is reduction of the matrix sizes in the 
associated Dirac operators, which leads to reduction of numerical costs in lattice QCD simulations. 
Thus they may be able to overcome the usual naive-fermion-based lattice fermions in lattice QCD 
[@]. The purpose of this work is reveal properties of staggered Wilson fermions in terms of the 
parity phase structure (Aoki phase) The Aoki phase for the staggered-Wilson was first studied 
in Ref. [||] and the present paper shows further investigation of this topic. The existence of the Aoki 
phase and the second-order phase boundary in Wilson-type lattice fermions indicates that one can 
apply them to lattice QCD simulations by tuning a mass parameter to take a chiral limit. Besides, 
the understanding of the parity-broken phase gives practical information for the application of its 
overlap and domain-wall versions. 

In this paper we elucidate the parity phase structure for staggered-Wilson fermions in the 
framework of the Gross-Neveu model and the hopping parameter expansion in the strong-coupling 
lattice QCD. We find the gap equations derived from the both theories show the pion condensate 
becomes nonzero in some range of the parameters and the pion becomes massless on the phase 
boundaries. It means the Aoki phase exists and the order of the phase transition is second-order. 
We also show we can take the chiral continuum limit in the Gross-Neveu model by tuning the mass 
and the gauge-coupling. These results on the staggered-Wilson fermion incidate we can obtain 
one- or two-flavor fermions by tuning the mass parameter and perform the lattice QCD simulation 
with these fermions as in the Wilson fermion. We note the results on the Gross-Neveu model is 
based on the work by some of the present authors [^, |6|] while the results on the strong-coupling 
lattice QCD are parts of a work in progress. 

2. Staggered Wilson fermions 

We begin with staggered-Wilson fermions in which the flavored-mass terms split the four 
degenerate tastes in a manner similar to the usual Wilson term. There are two possible types of the 
flavored-mass terms for staggered fermions as 

Mf = e £ f]ir]2mr]4CiC 2 C 3 C 4 = (1® /s) + 0(a), (2.1) 

sym 

M f = L -^v%T] v (C M C v + C v C M ) = (10 £ <v) + (2-2) 

n>v 2v J n>v 

whereC^ = (t^ = (-l) Xl+ - +x ^5 x , y , (e)^ = (-l)* + " +JC *& ,, (e^ = (-1)^4,; 

with (Vjjjry = U^ }X 5y yX+ ^. In the right hand sides we use the spin-taste representation as 1 75. We 
refer to M 1 as the Adams-type and the Hoelbling-type. The former splits the 4 tastes into two 
with positive(ra = + 1 ) and the other two with negative(ra = — 1 ) mass while the latter split them into 
one with positive(m = +2), two with zero(m = 0) and the other one with negative mass(»i = —2). 
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Now we introduce the Wilson parameter r = r8 X) y and shift the mass for the actions as with Wilson 
fermions. Then the Adams-type staggered-Wilson fermion action is given by 

Sa = £^[%0 M + r(l +M [ p) +M] xyX y, (2-3) 

xy 

with = \{Vn — V-^). Here M stands for the usual taste-singlet mass (M = M8 X y). The 
Hoelbling-type staggered-Wilson fermion action is given by 

Sh = £ Xx [%£>m + r(2 + m} 2) ) + M\ xyXy . (2.4) 

xy 

In the QCD simulation we will tune the mass parameter M to take a chiral limit. For some negative 
values of the mass parameter: — 1 < M < for Adams-type and — 2 < M < for Hoelbling-type 
with r = 1, we obtain two-flavor and one-flavor overlap fermions respectively by using the overlap 
formula. 

The potential problem in lattice QCD with these fermions is the breaking of some discrete 
symmetries as the shift symmetry caused by the flavored-mass terms [[I], Q]. There has not yet 
been a consensus on whether it does harm to lattice QCD with staggered-Wilson fermions. We 
can answer this question partly by studying the Aoki phase since a clear symptom is expected to 
appear in the phase structure if the symmetry breaking ruins the essential properties of QCD. In 
the following sections we will find the Aoki phase structure in the staggered-Wilson fermion is 
qualitatively similar to the original Wilson one and there is no disease. 

3. Gross-Neveu model 

We first investigate the parity phase diagram for staggered-Wilson fermions by using the d = 2 
Gross-Neveu model as a toy model of QCD. To study the pion condensate we generalize the usual 
staggered Gross-Neveu model to the one with the 75 -type 4-point interaction, which is given by 

; L^^ter+M -Xn-n)+Y,Xn(M+r(l +M f ))Xn 



2 ■ 



(3.1) 

A A J 

where the two-dimensional coordinate is defined as n = 2jV +A with sublattices A = (A\ ,A2){A\ 2 = 
0, 1). In this model Xn is a A^-component one-spinor (Xn)j(j =1)2, ...,A0 where XX = T!j=i XjXj- 
(_1 ^i+a 2 corresponds to I^ = 75 (g> 75 in the spinor-taste expression while = (— 
corresponds to 7^. In this dimension the Adams-type and Hoelbling-type flavored-mass terms co- 
incide and there is only one type Mf = r^r^ ~ 1 ® 75 + 0(a) with Ts = — iTJi^LsymCiCa. This 
mass term assigns the positive mass (m = +1) to one taste and the negative mass (m = — 1) to the 
other. With bosonic auxiliary fields o>-, n 'js leading to a-meson and 7r-meson fields, the action is 
rewritten as 

S = l^MXn+H -Xn-n)+Y*XnMfXn 



2 ■ 

n,n 



+ ^L((^-l-^) 2 + ^)+ £ X2^ + A(0^ + i(-\) M+Al K^) X 2^ + A, (3-2) 

Z 8 j¥ ,yVA 
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Figure 1: Aoki phase structure for the staggered- Wislon fermion in the Gross-Neveu model. A stands for a 
parity symmetric phase and B for Aoki phase. 



where we take r = 1 as the Wilson parameter. After integrating the fermion field, the partition 
function and the effective action with these auxiliary fields(meson fields) are given by 

Z = \ @a^®Kjre- m ^\ 5 eff =-^£((a^-l-M) 2 + ^)-TrlogD, (3.3) 

J A? jy 

withD„, m = (o> + i{-l) A,+A2 Ksr)8n,m + ! T($n+Li,m- 8n-n,m) + (.M f ) n!m . In the large N limit the 
partition function is given by the saddle point of the action as Z = e~ s ^ Cti ' no ' with the translation- 
invariant solutions ao, 7Co satisfying the saddle-point equations Ss < :if ^ 7Z °') = s ^ ff (ao,^a) _ q After 
some calculation process to derive the fermion determinant [||] we obtain the concrete forms of the 
saddle-point equations in the momentum space 



Cb-l-M / dk z (7 ((7 2 + ^o + s 2 ) -c\c\o ( ) 



(2n) 2 ((ob + cic 2 ) 2 + k 2 + s 2 ){(o - cic 2 ) 2 + n 2 + s 2 ) 



(3.4) 



TCo _ . f Jk 2 7r (g ( 2 + 7rg + s 2 ) + c\c\hq 

g 2 ~ J (2n) 2 ((a + c 1 c 2 ) 2 + n 2 + s 2 )((a -c l c 2 ) 2 + n 2 + s 2 )' 1 ' ' 

with = cosfc^/2 and = sinku/2. Now what we are interested in is the parity phase diagram 
in this theory. The parity phase boundary M c {g 2 ) is derived by imposing 71q = in (3_ L 4)( 3^5) after 
the overall 7To being removed in the second one. Then the gap equations are given by 

1+M, f dk 2 2c 2 c 2 o Q n » 

(3.6) 



g 2 J (2n) 2 {{o a + c l c 2 ) 2 + % 2 + s 2 ){{oo-c l c 2 ) 2 + K 2 + s 2 y 
1 _ „ f dk 2 a 2 + s 2 + c 2 c 2 



(27T) 2 ((ffo + CiC 2 ) 2 + K 2 + s 2 ){{o - CiC 2 ) 2 + 7T 2 + s 



By removing ao in these equations, we derive the phase boundary M c (g 2 ). The result is shown 
in Fig. [l[ It indicates the parity phase structure in the staggered-Wilson fermion is qualitatively 
similar to the usual Wilson case ^ reflecting the mass splitting of tastes given by the flavored 



mass. We also check the pion mass becomes zero on the second order phase boundary as m K oc 
V -jt s t\m=m c = where S e ff = VS e ff with V being the volume. 

O-7t 

We next consider the chiral and continuum limit of the staggered-Wilson Gross-Neveu models. 



The strategy is to expand the fermion determinant in the effective potential in Eq. (3.3) with respect 
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to the lattice spacing a. After some calculations (See details in Ref. [Q]) we obtain the effective 
potential remaining in the limit a — > 0, 

/M+l/a 2 \ / 1 ~ 1 7 \ , 



2^ ?! 



+ ( TT " Co + 2C 2 + - log 4a 2 ) a 2 + -(a 2 + tt ( 2 ) log ^°±^o . (3 . 8 ) 

with the three numbers as C = 1.177, C\ = -0.896 and C 2 = 0.404. Here taking the chiral limit 
means restoring the rotational symmetry in Ob and Hq by tuning the parameters. In this model we 
need introduce two independent coupling constants g 2 and g\ to restore the symmetry although 
the necessity of two couplings is just a model artifact. The tuned point for the chiral limit without 
0(a) corrections is 

M=-— G-l s 2 = — (3 9) 



To take the continuum limit we introduce the A-parameter as 2c?A = exp ^ Co — KC2 
the coupling renormalization for the chiral and continuum limit is given by 



Then 



^Co-^ + Ilog^), ^.ft + Ih,^). (3. 10 ) 

where we keep A finite when taking the continuum limit a — > 0. Finally the renormalized effective 
potential in the chiral and continuum limit is given by 

5eff = -(a 2 + 7r 2 )log^4? i ' ( 3 - n ) 

This wine-bottle potential yields the spontaneous breaking of the rotational symmetry. We have 
shown that the chirally-symmetric continuum limit can be taken by fine-tuning a mass parameter 
and two coupling constants in the staggered- Wilson Gross-Neveu model. Considering that the 
necessity of the two coupling constants is just a model artifact, this result indicates we can take a 
chiral limit by tuning only the mass parameter as in the Wilson fermion. Our results on the chiral 
and continuum limit for staggered-Wilson are qualitatively the same as the Wilson case 

4. Strong-coupling QCD 

In this section we investigate the Aoki phase structure in lattice QCD with staggered-Wilson 
fermions in the framework of the hopping parameter expansion (HPE) in the strong-coupling 
regime. We can detect a symptom of symmetry breaking from this analysis although we can- 
not know details of the true vacuum from this. For simplicity we concentrate on the Hoelbling-type 
lattice fermion here, but we can also make the same analysis in a parallel way for the Adams- 



type fermion. To perform the HPE for the Hoelbling-type fermion, we rewrite the action ( [2.4] ) by 
redefining X -> V^KX with K= 1/[2(M + 2r)], 

S = Y^XxXx + 2KY,Xx{rin D nhXy + IKr^UM^xyXy (4-1) 

x x,y x,y 
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Figure 2: (Left):Feynman rules for the HPE. (Center):one point function. (Right):two point function. 

In Fig. U(Left) we write down the Feynman rules in the HPE for this fermion. In this paper we 
perform the hopping parameter expansion up to 0(K ), which works for a small K. We derive 
chiral and pion condensates from the one-point function of the meson operator {jM x = XxXx) i n tne 
mean-field approximation. The equation for the one-point function up to 0(K 3 ) is obtained in a 
self-consistent way as shown in Fig. ^(Center), 

-I, = (^ x ) = (^)o + 2K 2 £ I, +A I, -2~ (Kr) 2 £ Z x Z x+fl+ v , (4.2) 

where we drop the link variable since we work in the strong-coupling limit. 0(K 3 ) diagrams are 
found to vanish due to the cancellation between diagrams. Here we solve this equation for the 
condensate £ within the mean-field approximation. For our purpose we assume L x = a x + ie x 7t x , 
where a x and 7i x correspond to the chiral and pion condensates. We substitute this form into Eq. 
( fk2l ) and obtain the self-consistent equation for o x and n x as 

-(o + ie x n) = -l + 2K 2 -4(o 2 + n 2 )-2~(Kr) 2 -4-3(o + is x x) 2 , (4.3) 

which yields — a = — 1 + \6K 2 K 2 and —in = — %K 2 -2io%. Here we have set r = 2a/2 for simplicity. 
We have two solutions depending on whether n = or n ^ 0: For % = we have a trivial solution 
a = 1. For 71 we have a non-trivial solution as 

CT = !6^ n = ± ^-{ X -^\ <4 ' 4) 

In this solution the pion condensate is non-zero and the ± signs indicate spontaneous parity break- 
ing. This parity-broken phase (Aoki phase) appears in the range of the hopping parameter or the 
mass parameter as | K |> 1/4 or equivalently — 4\/2 — 2 < M < —4\/2 + 2. We expect that the 
expansion up to 0(K 3 ) works to give a meaningful result at least around the critical parameter 
\K C \ 3 = 1/64 < 1. 

We next discuss the two-point function of the meson operator J^(0,x) = ^o^ x up to the 
same order as the one-point function. From Fig. ^|(Right) we derive the following equation for two 
point function. (0(K 3 ) diagrams vanish again.) 

^(o,x) = (xSx§x b x Xx) = -^ c +K 2 Y j {r (l x a (l x b x x b x ) 

+ { 2Kri ¥7l) LjxUy4+^ b x Xx)- (4-5) 
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Then the self-consistent equation for is given in the momentum space as 

2 

y(p). (4.6) 



K 2 V (e-^ + e ip " ) + ( 2Kr—^] V V e '^ ±Pv 



y(p)=-N c + 

We finally obtain the meson propagator as 

y(j>)=N c ' 



-2K 2 y cos p u +4 ( 2Kr = | Vcospncospv 



(4.7) 



The pole of 5?{p) gives the meson mass. Remembering 75 in the staggered fermion is given by 
e x = (— i)*i+-+*4 and the pion operator is given by % x = % x i£ x %x, it is obvious that the momentum 
of the pion should be measured from the shifted origin p = (n, %, %, 7i). Thus we set p = (im K a + 

\-\6K 2 



71, n, 71, 71) for 1/ y{p) = in ( |4.7[ ), which gives the pion mass m n as cosh(m^a) = 1 + 



6K 2 ' 

In this result the pion mass becomes tachyonic in the range | K |> 1/4. It indicates there occurs 
a phase transition between parity-symmetric and broken phases at |.ST C | = 1/4, which is consistent 



with the result of condensates in Eq. ( |4.4[ ). 

In this section we investigate parity phase structure for staggered-Wilson fermions by using 
the hopping parameter expansion up to 0(K 3 ) in the strong-coupling regime. Although we cannot 
give a strong argument just from this approximate calculation, we note that our result on the phase 
structure is qualitatively consistent with that of HPE for Wilson fermions which implies exis- 
tence of the Aoki phase in staggered-Wilson fermions. The full-order calculation in HPE and the 
effective potential analysis in progress [^] will give us more conclusive evidence for this topic. 

5. Summary 

In this paper we study the Gross-Neveu model and the strong-coupling lattice QCD with stag- 
gered Wilson fermions with emphasis on the Aoki phase structure. We have shown the parity 
broken phase and the second order phase boundary exist in the staggered-Wilson fermions as with 
the Wilson fermion. Our results indicate that we can apply the staggered Wilson fermions to lattice 
QCD simulations by mass parameter tuning. These results also indirectly suggest the applicability 
of the staggered overlap and staggered domain-wall fermions to lattice QCD. We note our results 
on the Aoki phase diagram exhibit no diseases due to a discrete symmetry breaking, which is 
consistent with the results in the lattice perturbation in [Q, W. 
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